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ABSTRACT 

We present a 2D study of the internal extinction on (sub)kiloparsec scales of a sample of local (z < 0.1) LIRGs (10) and ULIRGs (7), 
based on near-infrared Paof, BrS, and Bry line ratios, obtained with VLT-SINFONI integral-field spectroscopy (IFS). The 2D extinction 
(Av) distributions of the objects, map regions of ~ 3 X 3 kpc (LIRGs) and ~ 12 x 12kpc (ULIRGs), with average angular resolutions 
(FWHM) of ~0.2 kpc and ~0.9 kpc, respectively. The individual Ay galaxy distributions indicate a very clumpy dust structure already 
on sub-kiloparsec scales, with values (per spaxel) ranging from Av~l to 20 mag in LIRGs, and from Av~2 to 15 mag in ULIRGs. As 
a class, the median values of the distributions are Av=5.3 mag and Av=6.5 mag for the LIRG and ULIRG subsamples, respectively. 
In ~70% of the objects, the extinction peaks at the nucleus with values ranging from Av~3 to 17 mag. Within each galaxy, the 
Ay radial profile shows a mild decrement in LIRGs within the inner 2 kpc radius, while the same radial variation is not detected in 
ULIRGs, likely because of the lower linear scale resolution of the observations at the distance of ULIRGs. We evaluated the eff'ects 
of the galaxy distance in the measurements of the extinction as a function of the linear scale (in kpc) of the spaxel (i.e. due to the 
limited angular resolution of the observations). If the distribution of the gas/dust and star-forming regions in local LIRGs (63 Mpc, 
40pc/spaxel on average) is the same for galaxies at greater distances, the observed median Ay values based on emission line ratios 
would be a factor ~ 0.8 lower at the average distance of our ULIRG sample (328 Mpc, 0.2 kpc/spaxel), and a factor ~0.67 for galaxies 
located at distances of more than 800 Mpc (0.4 kpc/spaxel). This distance eff'ect would have implications for deriving the intrinsic 
extinction in high-z star-forming galaxies and for subsequent properties such as star formation rate, star formation surface density, 
and KS- law, based on Ha line fluxes. If local LIRGs are analogues of the main-sequence (MS) star-forming galaxies at cosmological 
distances, the extinction values (Ay) derived from the observed emission lines in these high-z sources would need to be increased by 
a factor 1.4 on average. 
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1. Introduction 

Since the first r esults obtained by the Infrared Astronomical 
Satellite (IRAS) ( [Soifer et"aL]|1984| ), there has been strong ef- 
fort to study the physical processes that power the luminous 
(LIRGs; 10i^Lo<LiR< IO^^l^) ^^^ ultraluminous (ULIRGs; 



% 



' 10^^Lo<LiR< 10^-^ Lo) infrared galaxy population (Sanders & 
|Mirabel|[Tf96l [Lonsdale et"aL]|2Q06j . The origin of the mid 



and near-infrared studies in LIRGs and ULIRGs dOenzel et al.| 
[T9981 [Scoville et al.||2Q00l [Alonso-Herrero et al.||20Q6| ) reveal 
that the distribution of the dust in these object is not uniform and 
that, though the dust tends to concentrate in the inner kilopar- 
secs with average visual extinction of Ay -3-5 mag in LIRGs and 
even higher in ULIRGs, the global distribution sh ows a patchy 
structure on kiloparsec and sub -kiloparsec scales ([Colina et al. 



and far-infrared emission (LiR[8-1000//m]) that dominates their 
bolometric luminosity is established as mainly due to massive 
starbursts with a small AGN contribution for LIRGs and with an 
increasing contribution in ULIRGs (e.g. 
Veilleux et al.||20Q9l [N^ini et al.|[JQTUl 



Goldaderet al. 1995, 



Alonso-Herrero et al., 



2Q12| and references therein). The radiation that originates in the 
starburst and/or the active galactic nucleus is then reprocessed 
by a surrounding dust component, and then re-emitted at long 
wavelengths in the form of a huge infrared emission. 

One of the main difl&culties in understanding the underlying 
power source of the LIRGs and ULIRGs is the high opacity of 
their nuclear regions. Previous optical ( [Garcia-Marin et al.|2009| ) 
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20001 [Garcia-Marin et al.|20Q6l|Bedregal et al.|2009| ) 

Besides the importance of knowing the 2D structure of the 
dust to understand the environment where the power source of 
the LIRGs and ULIRGs is embedded, dust plays a key role in the 
derivation of other physical and structural paramet ers of these 
objects, such as the derived star formation rate ( Garcia-Marin 
|et al.||2009| ), the eflTective radius fArribas et al.||2012| ), and as a 
consequence, the dynamical masses. 

Understanding the distribution and eflfect of dust in star- 
forming galaxies is also important for correctly interpreting or 
comparing diflferent tracers of star formation during the history 
of the Universe. This is in turn relevant when comparing lo- 
cal and high-z star- forming galaxy populations, which are often 
observed using diflferent tracers and/or resolutions. The distribu- 
tion of dust can in principle be studied in detail in local U/LIRGs 
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Table 1: The SINFONI sample 



IDl 


ID2 


z 


Dl 


Scale 


logLiR 


Common 


IRAS 




(Mpc) (pc/arcsec) 


(Lo) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


IRASF 12115-4656 


IRAS 12115-4657 


0.018489 


84.4 


394 


11.10 


IC5179 


IRAS 22132-3705 


0.011415 


45.6 


216 


11.12 


NGC 2369 


IRAS 07160-6215 


0.010807 


48.6 


230 


11.17 


NGC5135 


IRAS 13229-2934 


0.013693 


63.5 


299 


11.33 


NGC 3110 


IRAS 10015-0614 


0.016858 


78.4 


367 


11.34 


NGC 7130 


IRAS 21453-3511 


0.016151 


66.3 


312 


11.34 


ESQ 320-G030 


IRAS 11506-3851 


0.010781 


51.1 


242 


11.35 


IRASF 17138-1017 


IRAS 17138-1017 


0.017335 


75.3 


353 


11.42 


IC 4687 


IRAS 18093-5744 


0.017345 


75.1 


352 


11.44 


NGC 3256 


IRAS 10257-4338 


0.009354 


44.6 


212 


11.74 


IRAS 23128-5919 


IRAS 23128-5919 


0.044601 


195 


869 


12.04 


IRAS 21130-4446 


IRAS 21130-4446 


0.092554 


421 


1712 


12.22 


IRAS 22491-1808 


IRAS 22491-1808 


0.077760 


347 


1453 


12.23 


IRAS 06206-6315 


IRAS 06206-6315 


0.092441 


425 


1726 


12.31 


IRAS 12112-^0305 


IRAS 12112-^0305 0.073317 


337 


1416 


12.38 


IRAS 14348-1447 


IRAS 14348-1447 


0.083000 


382 


1575 


12.41 


IRAS 17208-0014 


IRAS 17208-0014 


0.042810 


189 


844 


12.43 



Notes. Col. (3): redshift from the NASA Extragalactic Database 
(NED). Cols. (4) and (5): Luminosity distance and scale from Ned 
Wright's Cosmology Calculator (Wright 2006 ) given ho = 0.70, Qm = 
0.7, Qm = 0.3. Col. (6): LiR( 8-1000//m) calculate d from the IRAS flux 
densities /i2> fis, /eo ^^^ /loo ([Sanders et al.|2003] ), using the expression 
given in |Sanders & Mirabel] ( |1996| ). 



with the advantage of the relatively high linear resolution and 
S/N. These studies can, therefore, help us interpret observations 
of analogous high-z star- forming galaxies, for which such a level 
of resolution, and S/N is not attainable with current instruments. 
The present work is part of a series presenting new H- and K- 
band SINFONI (Spectrograph for INtegral Field Observations in 
the Near Infrared, |Eisenhauer et al.|2QQ3| ) seeing-limited obser- 



vations of a s ample of local LIRGs and ULIRGs. |Piqueras Lopez 
et al.| ( [2012| ) (hereafter Paper I) presented the atlas of the sam- 



ple, the data reduction, and a brief analysis and discussion of the 
morphology of the gas emission and kinematics. In this second 
paper in the series, we focus on the study of the 2D distribu- 
tion of the dust derived using the Bry/Br^ and Paa/Bry ratios for 
LIRGs and ULIRGs, respectively, whereas in Piqueras Lopez et 
al. 2013 (Paper III, in preparation), we will apply the results for 
the 2D dust structure to study both the overall star formation rate 
(SFR) and the kpc structure of the SFR surface density (Ssfr) of 
the galaxies of the sample. 

The paper is organized as follows. In Sections [2] and [3] 
we briefly describe the sample, observations, and data reduc- 
tion process, which are detailed in Paper I. The procedures 
for obtaining the emission and Ay maps are described in Sec- 
tion |4j and the results and analysis of the Ay maps and distri- 
butions are presented in Section [5] Finally, Section [6] includes a 
brief summary of the paper. Throughout this work we consider 
Ho =70kms-^ Mpc'^ Da = 0.70, Qm = 0.30. 



2. The sample 

The SINFONI sample is a subsample of a larger set of local 
LIRGs and ULIRGs described in |Arribas eraL] ([2008) that was 
observed with diflTerent IFS facilities. These facilities include 
optical and infrared IFS instruments in both hemispheres, such 
as INTEGRAL-fWYFFOS ( |Arribas et aLl [T99g) at the 4.2 m 
WilHam Hers chel Telescope VLT- VIMOS {Visible MultiObject 
Spectrograph, \LeFhYrQ et al.|2003]), PMA S {Potsdam Multi Aper- 
ture Spectrophotometer, |Roth et al.||20Q5i ), and SINFONI. The 



large sample (-70 sources) covers the whole range of LIRGs and 
ULIRGs IR luminosities and the diflTerent morphological classes 
observed in these objects. 

The sample used for the present study comprises a total of 17 
objects, 10 LIRGs, and 7 ULIRGs covering the luminosity range 
log(LiR/L0)= 11.10 - 12.43 (see Table[T]). It was selected to be 
representative of the diflferent morphological types of LIRGs and 
ULIRGs, from isolated galaxies to strongly interacting systems 
and mergers. The mean redshifts of the LIRG and ULIRG sub- 
samples are zlirgs = 0.014 and zulirgs = 0.072 and the mean 
luminosities are log(LiR/Lo)= 11.33 and log(LiR/L0)= 12.29, 
respectively. More details on the sample can be found in Paper 
I. 

3. Observations and data reduction 

The data were obtained in service mode between April 2006 and 
July 2008 using SINFONI on the VLT (periods 77B, 78B, and 
8 IB). The sample was observed in the K band (1. 95-2.45 yum) 
with a plate scale of OV 1 25 xOV250 pixel" ^ that results in an FoV 
of 8"x8" by a 2D 64x64 spaxel frame. The spectral resolution is 
R~4000, and the full-width-at-half-maximum (FWHM) as mea- 
sured from the OH sky line at 2.190yum is ~6.0 A with a disper- 
sion of 2.45 A per pixel. The observations have a typical resolu- 
tion of -0.63 arcsec (FWHM, seeing-limited) that corresponds, 
on average, to -0.2 kpc and -0.9 kpc for LIRGs and ULIRGs, 
respectively. 

Owing to the limited FoV, the data sample typically ~ 3 x 
3 kpc for the LIRGs and ~ 12 x 12 kpc for the ULIRGs subsam- 
ple. Some of the more extended systems were then observed in 
diflferent pointings to cover regions of interest like secondary nu- 
clei or star-forming complexes. For a detailed description of the 
observations, pointings, and integration times, see Paper I. 

The reduction process was performed using the ESO pipeline 
ESOREX (version 2.0.5) and our own IDL routines for the flux 
calibration. All the individual frames were corrected from dark 
subtraction, flat fielding, detector linearity, geometrical distor- 
tion, wavelength calibration and sky- subtraction. After this pro- 
cess, the cubes of those objects with several pointings were com- 
bined to build a final mosaic. 

The maps of the Paa, Bry, and Br^ lines were constructed by 
fitting a Gaussian profile on a spaxel-by- spaxel basis. We have 
developed our own routines, based on the IDL routine MPFIT 
( [Markwar dt 2009), to perform the fitting of the cubes in an auto- 
mated fashion. As described in Paper I, the data were binned us- 
ing the Voronoi method developed by Cappellari & Copin|p003| ) 
in order to maximise the S/N over the entire FoV. Each map was 
binned independently, since the S/N depends on the wavelength, 
as well as the spatial distribution of the emission, to achieve a 
minimum S/N on average in the whole FoV. This minimum S/N 
varies from object to object, and is typically between 15 and 25 
for the brightest line (Bry and Paof in LIRGs and ULIRGs, re- 
spectively), and between 8 and 10 for the weakest (Br^ and Bry 
for LIRGs and ULIRGs, respectively). Further details on the flux 
calibration and the individual S/N thresholds used in the Voronoi 
binning can be found in the Paper I of these series. 

4. Data analysis 

The 2D extinction / dust structure was derived using the Bry/Br^ 
and Paof/Bry line ratios for LIRGs and ULIRGs respectively. Al- 
though the Br^ line is detected in most of the ULIRGs, its S/N is 
not high enough to map the emission and, in most of the cases, it 
is not suflftcient to perform an integrated analysis of the emission. 
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Fig. 3: Surface brightness distributions of the individual spaxels of LIRGs and ULIRGs. For clarity, we only plot a random 
distribution of the 20% and 50% of the points for LIRGs and ULIRGs, respectively. The colour code indicates the Ay value of each 
spaxel, whereas the points with Av< are outlined with a black contour line. The increase in the extinction towards low surface 
brightness values is mainly due to the high uncertainties of the flux measurements, in particular in ULIRGs, where the highest Ay 
values correspond to the spaxels with the lower Bry surface brightness. 



As mentioned before, the maps of the difl'erent lines were 
constructed by fitting a single Gaussian profile on a spaxel-by- 
spaxel basis (see Fig. [T] and |2]). Based on the emission maps, 
we obtained the extin ction in magnitudes (A y) following the 
procedure outlined in [Bedregal et al.| ( |2009| ). We compared 
the theoretical ratio between the two lines (Br7/Br^= 1.52 and 
Fsia/Bry= 12.07 at T= 10, 000 K and ne = 10^ cm-\ case B; 



Qsterbrock||1989| ) with the measurements for each spaxel. The 
extinction in magnitudes could be expressed in the form 



■A, 



-2.5- log 



(F^^F^Jq 



(Fii/F^2)T 



(1) 



where F^.^o and F^. j are the observed and theoretical fluxes for 
a line centred at Ai. We made use of the extinction law described 
in Calzet ti et al.| ( |2000 ) to express Equation [T] in terms of the 
visual extinction Ay (Aery = 0.096 Ay, Abvs = 0.1 32 Ay and 

Apa, = 0.145 Ay). 

Since the individual values of Ay are sensitive to the S/N of 
the weakest line (Br^ and Bry for LIRGs and ULIRGs respec- 
tively), we have only considered those spaxels where the weakest 
line has been detected above an S/N threshold of four to obtain 
reliable Ay. This efl'ect is very significative in the case of the 
BryfBrS ratio since the Br^ line lies close to the blue limit of 
the SINFONI K-band. As discussed in Paper I, this wavelength 
region is strongly aff'ected by noise due to the sky emission, and 
the atmospheric transmission also decreases. This translates into 
a more complex local continuum determination, making the line 
fitting more uncertain. An excess in the continuum level esti- 
mation would decrease the line flux and, therefore, increase the 
extinction (see expression above). Although the PacK line also 
lies in this region of the spectra, this efl'ect is not so relevant, 
given the strength of the line (more than xlO the Bry emission), 
and it is in the numerator of Equation [T] 

The Icr uncertainties of the individual Ay values vary 
typically from 10-20% in central regions with high S/N, 
up to 70-80% in external areas of low surface brightness 



(S < 10"^^ erg s"^ cm"^ kpc"^), with a median value of 30-35%. 
Owing to the larger individual errors in the low surface bright- 
ness areas, we observed an artificial increase in the Ay mea- 
surements, especially in the ULIRGs. This systematic eff'ect is 
observed in Fig. [3] where the highest extinction values are mea- 
sured in those spaxels with the lower Bry surface brightness. The 
Ay uncertainties are obtained from the corresponding errors in 
the fluxes of the two lines, which as described in Paper I, are 
estimated using Monte Carlo simulations. The advantage of this 
kind of error estimation is that uncertainties are directly mea- 
sured from the spectra, so they not only include the efl'ect from 
photon noise but also take uncertainties due to an improper con- 
tinuum determination or line fitting into account. 

5. Results and discussion 

5.1. Two-dimensional extinction structure in LIRGs and 
ULIRGs 

Our 2D extinction maps cover areas from ~2kpc to ~12kpc 
for LIRGs and ULIRGs, respectively (see Figs. [T] and [2]). The 
seeing-limited observations provide a linear resolution equiva- 
lent to a physical scale resolution of -0.2 kpc and -0.9 kpc for 
each subsample. On these scales, the extinction maps show that 
dust is not uniformly distributed, revealing a clumpy structure 
with almost transparent areas (Ay < 1 mag) and regions where 
the visual extinction is higher than ten magnitudes. 

In LIRGs, the extinction maps show a very irregular and 
clumpy structure on scales of -200-300 pc, already observed 
from near-IR continuum maps ( [Scoville et al.|200"0| ). The higher 
Ay values are usually associated with the nuclear regions of the 
objects, although obscured extranuclear regions are also com- 
mon. The star-forming regions with high Bry surface brightness, 
found along the dynamical structures like arms of rings, are typ- 
ically low-extinction regions. 

As shown in the radial profiles of NGC 3110 or IRASF 
12115-4656 (Figs, lb and[Te]), the extinction increases inwards 
up to Ay- 15-25 mag. This behaviour is also observed in ESO 
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Table 2: Integrated properties and statistics of the Ay distributions 



Object 
(1) 



Av,ni 



(2) 



Av,nuclear Av,Rgjf 



(3) 



(5) 



Av,FoV 

(6) 



(7) 



(8) 



AvCPs) 
(9) 



Av (P95) 
(10) 



-^^ spaxel 

(11) 



IRASF121 15-4656 








7.7 + 0.1 


5.2 + 0.1 


4.9 


6.3 + 8.1 


-3.65 


21.5 


1533 


IC5179 


10.4 + 0.4 






4.2 + 0.1* 


4.2 + 0.1 


4.4 


5.2 + 6.3 


-3.00 


17.3 


2099 


NGC2369 


16.8 + 0.2 






17.3 + 0.1 


15.0 + 0.1 


15.5 


13.5 + 6.6 


3.75 


26.0 


603 


NGC5135 


2.8 + 0.2 






3.8 + 0.1 


3.7 + 0.1 


4.2 


4.7 + 4.5 


-1.72 


11.7 


1100 


NGC3110 








7.4 + 0.1* 


7.4 + 0.1 


7.9 


7.9 + 5.8 


-3.02 


16.6 


400 


NGC7130 


12.3 + 0.3 






11.8 + 0.2 


8.7 + 0.1 


6.2 


5.9 + 6.5 


-2.45 


18.6 


687 


ESO320-G030 








7.8 + 0.1 


7.0 + 0.1 


7.8 


8.0 + 7.1 


-2.82 


20.1 


1383 


IRASF17138-1017 


15.2 + 0.4 






11.1 +0.1 


7.0 + 0.1 


7.6 


6.1 + 5.6 


-0.91 


17.2 


806 


IC4687 


10.6 + 0.1 






3.9 + 0.1 


3.7 + 0.1 


4.2 


4.5 + 5.1 


-1.74 


13.1 


3401 


NGC3256^ 




12.2 + 0.2 


5.0 + 0.1 


4.1 +0.1 


5.0 


6.6 + 6.1 


-2.73 


16.8 


4522 


IRAS23 128-5919 


8.7 ± 0.3 


6.7 ± 0.2 


7.1 +0.1 


6.8 + 0.1 


6.2 


7.0 + 4.7 


-1.60 


13.4 


1182 


IRAS21 130-4446 


4.3 + 0.4 


3.2 + 0.2 


4.7 + 0.3 


4.4 + 0.1 


5.4 


5.0 + 5.9 


-1.80 


15.5 


214 


IRAS22491-1808 


4.1 ±0.3 


4.5 + 0.2 


5.5 + 0.1 


5.3 + 0.1 


6.4 


6.5 + 5.1 


-1.79 


16.0 


299 


IRAS06206-6315 


7.5 + 0.3 


11.3 + 0.4 


7.6 + 0.2 


8.5 + 0.2 


10.7 


9.5 + 7.7 


-3.05 


22.7 


114 


IRAS 121 12-^0305 


8.9 ± 0.3 


8.4 + 0.2 


9.2 + 0.2 


8.5 + 0.1 


8.2 


9.0 + 6.2 


-0.52 


20.9 


466 


IRAS 14348- 1447 


5.5 ± 0.2 


8.1+0.3 


6.2 + 0.2 


7.1 +0.1 


8.6 


7.9 + 6.9 


-1.04 


20.6 


301 


IRAS17208-0014 


8.0 + 0.2 






6.8 + 0.1 


5.8 + 0.1 


4.8 


6.3 + 4.0 


-1.61 


10.8 


485 



Notes. Cols. (2) and (3): Nuclear extinction measured at the main ( 2) and secondary (3) nucleus within an aperture radius of 0V63. Cols. (4) 
and (5): Measurements of the Ay within the Ha effective radius from |Arribas et al.]p012| ), and the integrated measurement over the whole FoV, 
respectively. Cols. (6) to (9): Median, weighted mean, and 5th and 95th percentiles of the Ay distributions. Col. (10): Number of valid spaxels in 
the Ay maps and distributions. * Since the main nucleus of NGC 3256 was not observed, we centred the aperture for Ay r^^ on the centre of the 
FoV, so the measurement might be inaccurate. *Due to the limited FoV of the observations, in these objects the Ha effective radius is greater than 
our FoV, and the Ay r^^, measurements are equivalent to Ay pov- 
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Fig. 4: Individual Ay distributions, on a spaxel-by- spaxel basis, 
of the galaxies of the sample ordered by increasing Ljr. The 
extremes of the distributions are the 5th and 95th percentiles (P5 
and P95). The boxes illustrate the interquartile range, whereas the 
horizontal blue and red lines correspond to the median and the 
weighted mean of the distribution, respectively. The measure- 
ments of the nuclear extinction are plotted as black diamonds. 
The total Ay distributions for LIRGs and ULIRGs are shown on 
the right-hand side of the plot. 



320-G030 (Fig. Id), where the emission is concentrated on a 
star-forming ring of -500-600 pc radius. Although the nucleus 
might also be obscured, the radial profile of this object is not as 
steep as in NGC 3110 or IRASF 12115-4656 and indicates that 
the lack of emission could also be due to the intrinsic distribution 
of the star- forming regions around the ring. 



The morphology of the Ay maps in the ULIRG subsample 
suggest a patchy, non-uniform distribution of the dust, typically 
on physical scales of > 1 kpc that correspond to our resolution 
limit. Owing to the higher linear resolution (i.e. kpc/spaxel) of 
the ULIRG subsample, the comparison with the LIRGs is not 
straightforward. As shown in Fig. |5j although our data sam- 
ples similar areas of ~l-2Reff, the dust structure is probed with 
significantly different spatial resolutions, owing to the factor x5 
in distance between both subsamples. This difference precludes 
our from resolving sub-kiloparsec structures in ULIRGs, such 
as the ones observed in the LIRG Ay maps, limiting our physi- 
cal resolution to ~1 kpc. As discussed in Sec.|5.4[ these differ- 



ences in the linear resolution between both subsamples not only 
shape the observed dust morphology in the more distant galaxies, 
but also might determine global measurements of the extinction, 
such as the median of the Ay distributions. In Sec. [53] we dis- 
cuss how this distance effect would have direct implications for 
the study of high-z galaxies. 

5.2. Ay distributions and radial profiles 

Figures [T] and [2] show the Ay distributions for each galaxy. Al- 
though it is clear that most of the spaxels with Ay<0 have no 
physical meaning individually, we have kept them in the distri- 
butions since they do have statistical relevance. If we remove 
them from the distributions, we introduce a bias toward the high 
Ay values, displacing the mean and median of the distributions 
artificially. On the other hand, due to the S/N threshold adopted, 
we assure that most of those spaxels with Ay<0 are compatible 
with Ay~0 within the uncertainties. 

The histograms show a wide variety of distributions, from 
narrow, peaked distributions, such as NGC 5135 or IC 4687 
(Fig. [T^and[Th]), concentrated towards low Ay values, to wide 
distributions such as NGC 3110 or IRAS 14348-1447 (Figs.flb 



and[2c]) that extend up to -30-35 mag. The median and weighted 
mean Ay values, together with the 5th and 95th percentiles of the 
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Fig. 5: Av distributions and radial profiles of the LIRGs (top) and ULIRGs (bottom) subsamples on a spaxel-by-spaxel basis. The 
median Ay values of each distribution are shown in the left panels. In the central and right panels, the radial profiles are plotted in 
terms of the radius in kpc (centre) and in units of the Ha eff'ective radius (Reff, right panel) extracted from Arribas et al., ( ,2012) . The 
red line represents the weighted mean of Ay and its error for diff'erent radial bins in steps of 1/30 of the total radial coverage. For 
the ULIRG subsample, we plotted the radial profile of the sample by considering each component of the systems separately (top) 
and the profile extracted by taking the brightest nucleus in the K-band continuum as the centre of the systems (bottom). 



distributions, are listed in Table[2] Figure[4]also compares the in- 
dividual distributions of each galaxy of the sample, ordered by 
increasing Ljr. As shown in the figure, most of the individual 
values, within the interquartile ranges, are concentrated between 
Ay~l and Ay~20 mag, and there is no clear evidence of any 
dependence with Ljr. 

In LIRGs, the visual extinction ranges between Ay~ 1 - 
20 mag, whereas in ULIRGs, the Ay values range between Ay~ 
2-15 mag. In LIRGs, these values of the visual extinction are 
very similar although slightly lower than previous results in the 
mid-infrared from Spitzer, Ay < 1 - 30 mag with a mean value 
of -11 mag ( [Pereira-Santaella et al.|2010]|Alonso-Herrero et al.| 
[2012). On the other hand, ULIRGs show slightly higher values 
than previo us measurements in the op tical, from Ay < 0.2 mag 
to ~ 9 mag ( |Garcia-Marin et aL]|2009| ), and significantly lower 



than mid-infrared measurements 



3ased on the silicate absorption 



feature at 9 .7 yum, from Ay~6mag up to Ay > 40 mag (Imanishi 
let al.|2007| ). 

These Ay values have to be considered as lower limits of 
the dust extinction, since the theoretical values of the ratios are 
based on the assumption that the gas is optically thin. However, 
it is well known that the central regions of these objects are dusty 
environments and that the observed recombination lines might 
have been originated at diff'erent depths within the star-forming 



regions. That also means that we would infer diff'erent extinc- 
tion values from diff'erent emission line ratios, since the lines are 
probing diff'erent optical depths. 

To compare the Ay distributions of LIRGs and ULIRGs, we 
combined all the available spaxels of each luminosity bin to ob- 
tain a typical Ay distribution. These distributions are shown in 
Fig. [4] and, in more detail, in Fig. |5] The median values of Ay 
for each subsample are similar, Aymed = 5.3 mag and Aymed = 
6.5 mag, and the number of individual spaxels is -16000 for the 
LIRG distribution and ~3000 for the ULIRG one. The shape 
of the distributions are, however, slightly diff'erent. Whereas the 
LIRG distribution seems to have more than 50% of the points 
concentrated within the range Ay^ 1 - 10 mag, the ULIRGs 
extend over a somewhat narrower range of Ay^ 3 - 10 mag 
and tend to reach higher Ay values on individual spaxels. On 
the other hand, the modes of the distributions are also diff'er- 
ent; whereas the LIRG distribution peaks at ~ 3-4 mag, the 
mode of the ULIRG distribution reaches up to ~ 7 - 8 mag. Al- 
though some of these diff'erences might be intrinsic, we discuss 
in Section |5.4| that the physical sampling of the maps plays an 
important role in the study of the extinction, owing to the patchy 
structure of the dust. 

We obtained the radial profiles of the extinction for every in- 
dividual object of the sample and characteristic profiles for both 
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LIRG and ULIRG subsamples. For this purpose, we adopted the 
same criterion as in Paper I to identify the central spaxel of each 
object with the brightest spaxel of the FoV in the K-band image 
(see Figs [T] and |2]). The only exception in NGC 3256, since the 
nucleus of this galaxy was not observed in the K-band (see Pa- 
per I for further details). For this object, we used the H-band 
continuum image of the nucleus to identify the central spaxel. 
For those systems with multiple nuclei (i.e. all the ULIRGs with 
the exception of IRAS 17208-0014), we also extracted the radial 
profiles of each component separately. 

As shown in Figs. [T] and [2j the profiles typically sample the 
inner ~ 2 kpc for the LIRGs and < 6 kpc for the ULIRGs, and 
most of them have an almost flat or negative slope. LIRGs 
show steeper negative slopes than ULIRGs, especially in the 
central 0.5-1 kpc. In LIRGs, these slopes are typically of 



Distance (Mpc) 
500 1000 



-2.4 mag kpc on average, versus 



-0.3 mag kpc in 



ULIRGs. This could be explained by the difl'erent sampling 
scales of both subsamples, since we cannot resolve the innermost 
regions of the ULIRGs with the resolution achieved for the LIRG 
subset. In those pre-coalescent systems with multiple nuclei, we 
found no systematic difl'erences between the radial profiles of 
both components. These profiles typically sample the innermost 
-2-3 kpc of each component of the system, and show no clear 
radial dependence of the extinction on these scales. 

We extracted the radial profiles of the total Ay distributions 
of LIRGs and ULIRGs separately. Since each object has been 
observed with diff'erent sampling and physical scales, we also 
obtained the profiles in units of the eff'ective radius Reff , using 
the values fro mjArribas et al.| ( |2012| ) obtained from Ha maps. As 
shown in Fig.[5| the LIRG profile is almost flat beyond r~l kpc 
or -0.5 Reff, with an average value of Av~5.3 mag. Within the 
central kiloparsec, the extinction increases up to ~ 10 mag. The 
ULIRG subsample shows a very uniform profile, with a median 
value of Ay -6.0 mag and only local deviations due to the pres- 
ence of double nuclei in some of the systems, or due to strong 
complexes of star formation at distances beyond 2-3 kpc radius 
(see Fig|2a|andl2b]for some examples). The radial profile shows 
small diff'erences when extracted for each component separately, 
with an almost flat slope over the inner 2-3 kpc radius (~Reff). 
The measurements beyond r ~6 kpc or r ~2 Reff are very uncer- 
tain, owing to the lack of available spaxels and to the low surface 
brightness of the Bry emission in these regions. As mentioned 
before, the extinction in ULIRGs derived using Eq. [T] is highly 
aff'ected by noise fluctuations of the Bry line. 

5.3. Nuclear and integrated Ay measurements 

We obtained integrated Ay measurements for difl'erent regions 
of interest in each object (Table [2]). The uncertainties of the pa- 
rameters are obtained by a Monte Carlo method of N = 1000 
simulations, and do not take the Icr uncertainties of ~5% into 
account in the absolute flux calibration (see |Piqueras Lopez et al.| 
[2012). We found that, in -70% of the objects, its nucleus cor- 
responds to a peak in the extinction, ranging from Ay- 3 mag up 
to Ay~17 mag. These values of the nuclear extinction are higher 
than each median and mean Ay values in ~ 57% of the objects. 
However, there are some galaxies, such as NGC 31 10 or IC 4687 



(Figs.[Tb]and[Th]), where the nucleus is completely obscured, and 
no measurements of the extinction are available. That the nu- 
cleus of the objects coincides with a peak of the extinction agrees 
with other studies in LIRGs and ULIRGs (see ' Alonso-Herreroj 
[eraLl[200 6 and Ga rcia-Mann et al.f2009j based on Ua/Ufi and 
Paa/Bry ratios, and with mid-infrared studies with Spitzer, based 
on the silicate absorption feature at 9.7 //m ( [Imanishi et al.|200T 
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Fig. 6: Evolution of the computed median Ay of the indi- 
vidual distributions of the LIRG subsample as a function of 
scale/distance. The Ay values are normalised to the rest- framed 
distribution. The value that corresponds to the Ay distribution 
of the whole LIRG subsample is plotted as black circles and a 
power-law fit to the data as a dashed black line. The mean dis- 
tance of the ULIRG subsample (328 Mpc) is plotted as a vertical 
dotted line for reference. 



[Pereira-Santaella et al. 112010] ). These authors found that either 
the highest extinctions coincide with the nucleus of the objects 
or the nuclear regions are local maxima in the Ay maps. 

For those objects with a double nucleus, we also presented 
integrated measurements of the extinction for the secondary one. 
Although almost all of the interacting systems in our sample are 
ULIRGs, NGC 3256 presents a well known, highly obscured 
nucleus ~ 5 arcsec to the south of the main one. Our measure- 
ments of the extinction of the southern nucleus of this object 
shows that it is one of the most extinguished regions in our sam- 
ple, with Ay ^12.2 mag, in good agreement with previous works 
CKotilainen et al.'1996"Alon so-Herrero et al.|2006||Diaz-Santosl 
,etal.2008 „Richetal..2011j . 

5.4. Dust dumpiness and tlie effect of tlie linear resolution. 



As mentioned in Section 5.1 the Ay maps reveal a patchy, 
clumpy structure of the dust at sub-kiloparsec and kiloparsec 
scales in LIRGs and ULIRGs, respectively. Given this non- 
uniform distribution, the measurements of the extinction at dif- 
ferent distances might be aff'ected by the physical scale of the ob- 
servations. To probe how the pixel scale might bias the Ay mea- 
surements, we have obtained the Ay distribution of our LIRG 
sample simulating diff'erent scales, hence diff'erent distances. 
This distance eff'ect because of the linear resolution might be 
even more relevant for high-z objects, where the structure is sam- 
pled on even larger scales of -1-2 kpc. 

To simulate the distribution at further distances, we degraded 
the individual Bry and Br^ maps to diff'erent scales, and obtained 
maps of poorer spatial resolution. In this process, we only con- 
sidered the same valid spaxels as in the original maps. Once the 
maps were degraded, we obtained the Ay distributions of each 
individual object as described in Section]?] which were finally 
combined in one single distribution. 
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Fig. 7: Different simulations of the LIRGs Ay distribution for increasing distance (from top to bottom and left to right). The median 
of each distribution is highlighted with a vertical dashed line. The distance and the median Ay of the simulated distribution is 
annotated in each panel. The top left panel corresponds to the average distance of the ULIRG subsample, and shows the observed 
Ay distribution of the ULIRGs and its median, in yellow. 



Since the FoV of our SINFONI observation is limited to ~ 
8" X 8", we could only sample the innermost ~ 3 x 3 kpc of the 
LIRGs (typically -Reff). If we translate these scales to a distance 
that is ten times larger, these ~ 3 kpc are equivalent to an angular 
distance of ~1", i.e. 1/8 of the FoV or ~ 8 spaxels. The lack 
of data from the external regions of the local objects keeps the 
extrapolation to larger distances from beeing straight-forward. 

Figure [6] shows the relation of the median of the individual 
Ay simulated distributions to the sampling scale/distance, for 
each galaxy in the LIRG subsample. Although each curve shows 
a different behaviour that depends on the particular gas and dust 
distribution, there is a general trend towards lower Ay values as 
distance increases. This decrement of the median value of the 
visual extinction is also observed when we consider the distri- 
bution of LIRGs as a class, adding all the individual spaxels of 
the LIRG subsample in a single distribution, as shown in the fig- 
ure. To parametrize this observed decrement of the median Ay 
of the LIRG distribution with the sampling scale/distance, we 
fitted the data points to a simple power-law model, and found 
that the best-fit model corresponds to 



Ay/Ayo - (D/Do)- 



-0.13 



(S/So) 



-0.14 



(2) 



where Ay and S are the median values of the Ay distribution 
and the physical scale at a distance D, respectively, and Ayo = 
5.27 mag and So = 290 pc arcsec"^ are the median values of the 
rest-frame Ay distribution and median physical scale at the mean 
distance of the LIRGs subset. Do = 63.3 Mpc, respectively. 

In figure [7] we show in detail the observed Ay distribution 
of the LIRG subsample and the simulated distributions for in- 
creasing distances. The different panels reveal that not only the 



median value of the distribution decreases when the galaxies are 
sampled on a larger physical scale, but also that the shape of 
the distribution is slightly different, becoming narrower than the 
rest-framed distribution, and more compact. The mode of the 
simulated distribution also changes with respect to the ULIRG 
observed distribution, and becomes lower than the median. As 
shown in Fig. [6j the difference between the rest-frame extinc- 
tion, Ayo, and the simulated Ay increases rapidly within the 
first D~700Mpc, and seems to reach an asymptotic value of 
Ay/Ayo -0.65 beyond that distance. 

The non-uniform distribution of the dust, even on small 
scales, makes that, at a given resolution unit, we map both ob- 
scured regions and areas where the interstellar medium is more 
transparent. This average is biased towards the brightest regions, 
which are those were the emission is less absorbed. The result 
is that, on average, the dust distribution is smoothed, becomes 
narrower, and the Ay values that we measure are biased towards 
the lowest values. 



As mentioned in Section [5^ we measured a median extinc- 
tion of Ay LIRGs = 5.3 and Ay ulirgs = 6.5 magnitudes for our 
LIRG and ULIRG subsamples, respectively. Since the average 
distance of each subsample is ~ 63 Mpc and -328 Mpc (see 
Paper I), this factor ~5 in distance and sampling scale (from 
-0.2 kpc to -0.9 kpc on average, respectively) translates to a de- 
crease in the measured Ay of ~L2mag, as shown in Fig. [7] (top 
right panel). If we assume that LIRGs and ULIRGs have similar 
structures and apply the same correction to the ULIRGs subsam- 
ple, we find that the average Ay in the ULIRGs, corrected from 
distance effects, reaches Ay~8.0mag, i.e. -2.8 mag higher than 
in LIRGs. 
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Besides this resolution/distance effect, the difference be- 
tween the observed Ay distribution in LIRGs and ULIRGs could 
also be interpreted in terms of intrinsic differences in the mor- 
phology between both classes. Although there seems to be a 
general trend for LIRGs as a class that suggest that the observed 
Ay decreases with the distance/resolution, Fig [6] also shows that 
the behaviour of each individual galaxy is very different, and 
depends strongly on the particular morphology of the dust distri- 
bution. These differences among the objects of the LIRGs subset 
suggest that the correction to the observed Ay could not be ac- 
curate for individual galaxies, and would only reflect a general 
trend for LIRGs as a class. 

5.5. Implication for extinction-corrected properties in high-z 
galaxies 

There are different observational proofs that suggest that sub- 
millimetre galaxies (SMGs) could be the high-z analogous of 
local (U)LIRG. The sub-millimetre and radio fluxes of this pop- 
ulation of galaxies indicate that their bolometric luminosities 
are comparable to local ULIRGs, whereas their mid-IR emis- 
sion seems to be similar to the observed in local LIRGs (Kovacs 
let al.|20Q6|[Takata et al.|20Q6l|Menendez-DelmesEeeral.|20Q9| ). 
Besides this, recent IFS -based observations of SMGs reveal that 
they present evidence of clumpy star formation on kiloparsec 
scales and similar star-formation rate surface densities (Ssfr) 



to the local counterparts ([Nesvadba et al^|2QQ7[ [Harrison et"aL 



2Q12[|Alaghband-Zadeh et al.|2Q12||Menendez-Delmestre et al. 



2013 1. It is well known that these high-z galaxies may be highly 
obscured and that a combination of intense dust-obscured star 
formation and dust-enshrouded AGN activity would be respon- 
sible for the high infrared luminosities of these objects ( |Blain| 
let al.|2QQ2) . HST-NICMOS and ACS observations reflect st ruc- 
tured dust obscuration i n these objects ( |Swinbank et al. 12010 ). In 
particular, |Takata et al.| ( |2006] ) find that the internal extinction in 
these objects are similar to the extinction in local ULIRGs, and 
measured a median extinction of Ay = 2.9 ± 0.5 mag in a sample 
of SMGs at z~ 1.0-3. 5 using the HJH^ flux ratio. 

Clumpy and dusty star-forming structures have also been 
identified at high redshifts in more common star forming galax- 
ies (i.e. the so called Main Sequence star forming galax- 
ies). These galaxies have (sub)kpc star-forming clumps mostly 
spread over galactocentric distances of few to several kpc 
( [Forster Schreiber et all2QTTal|2011b||Genzel et al.|2011| ), and 
with internal nebular extinctions of 2 to 4 magnitudes, assuming 
AyH, = Ay Stellar /0.44 ( |Forster Schreiber et al.|2009l|Wuyts et al. 
[Mfg. 

In the previous subsection we showed that the distance/linear 
scale may play an important role in deriving the internal extinc- 
tion properties of LIRGs and ULIRGs, by comparing both pop- 
ulations of galaxies locally, albeit with a difference of a factor 
x5 in distance between both subsamples. According to the sim- 
ulations, the median extinction measured in a given galaxy de- 
creases when placed at increasing distances. This effect is due 
to the fixed angular resolution of the IFS data that translates into 
a larger physical scale per spatial resolution element (spaxel) as 
the distance to the galaxy increases, and is particularly important 
when the physical scales sampled by the spaxel are much larger 
than that of the intrinsic clumpy structure of the dust distribution 
and star- forming regions. 

It is clear that if the intrinsic star- forming structure of high- 
z galaxies is in general similar to that of our local LIRGs (i.e. 
mostly disks) and ULIRGs (i.e. mostly interacting), the smear- 
ing effect mentioned above would have a direct impact in the 



derivation of their 2D internal extinction values, and of all rele- 
vant subsequent extinction corrected properties such as star for- 
mation surface densities, KS-laws, and overall star formation 
rates. This would certainly be the case not only with seeing- 
limited near-IR IFS, heavily undersampling the galaxies at red- 
shifts of 1 to 3, with each spaxel corresponding to about 1.5-2 
kpc, but even when using AO assisted IFS where the spaxel (50 
to 100 mas) translates to about 0.4 to 0.8 kpc. Thus, if the results 
given by Eq. |2] (see also Fig. [6]) are applied to high-z galaxies, 
the extinction values derived directly from the observed optical 
emission line ratios would require to be increased by a factor 
1.4 on average. This correction would correspond to an addi- 
tional increase in the H^, flux by a factor ~3 and hence, in the 
extinction-corrected SFR and Ssfr. 

Finally, it is worth noticing that this distance effect could 
also have a minor wavelength dependency. The method pre- 
sented here to derive the Ay values is based on specific emission 
line ratios in the near-infrared. Since these lines originate in re- 
gions of higher optical depths than the optical emission lines, 
and the stellar continuum measured at optical wavelengths, dif- 
ferences in the amount and/or distance dependency could appear. 
It would be worth exploring these effects with a suitable set of 
data, in particular if, as in many high-z studies, the Ay correc- 
tions applied to the observed Ha flux is obtained indirectly using 
the standard Calzetti recipe Ah^ = 7.4 E(B-V), where E(B-V)gas 
= 0.44E(B-V),,^,, (Calzetti et al.|2000|[20QT] ) 

6. Summary 

In this paper, we presented a detailed 2D study of the extinction 
structure of a representative local sample of 10 LIRGs and 7 
ULIRGs, based on VLT-SINFONI IFS K-band observations. We 
sample the central ~ 3x3 kpc for LIRGs and the ~ 12x 12 kpc for 
ULIRGs, with average linear resolutions (FWHM) of ~0.2kpc 
and -0.9 kpc, respectively. The extinction maps are based on 
measurements of the Bry/BrS line ratio for LIRGs and of the 
Paof/Bry line ratio for ULIRGs. 

In agreement with previous studies, we found that the dis- 
tribution of the dust in these galaxies presents a patchy struc- 
ture on sub-kiloparsec scales, with regions almost transparent 
with Ay~0 to heavily obscured areas with Ay values up to 
-20-30 mag. In most of the objects in the sample (-70%), the 
nucleus of the galaxy coincides with the peak in the extinction 
maps, with values that range from Ay~ 3 mag up to Ay ~ 17 mag. 

We obtained the Ay distribution of the individual galaxies on 
a spaxel-by-spaxel basis (see Fig.|4]). The individual Ay distri- 
butions show a wide range of values with most of them spread 
between Ay~l and Ay~20 mag, with no clear evidence of any 
dependence with Ljr. However, as a class (see Fig. [5]), ULIRGs 
show Ay values (median of 6.5 mag, mode of -7-8 mag) higher 
than those for LIRGs (median of 5.3 mag, mode of ~3-4mag). 
The Ay distribution in LIRGs shows a mild decrease as a func- 
tion of galactocentric distances of up to 1 kpc and a flattening at 
larger distances (2-3 kpc). No similar behaviour is detected in 
ULIRGs, most likely owing to the lower linear resolution of the 
observations. 

To study the effect of the spatial sampling (i.e. physical 
scale per spaxel) in the derived extinction values at increasing 
galaxy distances, the individual Bry and Br^ maps of our sub- 
sample of local LIRGs (at an average distance of 63 Mpc) have 
been artificially smeared. These simulations have shown that 
the spatial resolution plays an important role in shaping the Ay 
distributions. The median value of the visual extinction mea- 
sured on the LIRG subsample decreases as a function of the 
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linear resolution/distance by a factor ~ 0.8 at the average dis- 
tance (328 Mpc, 0.2 kpc/spaxel) of our ULIRG sample, and up 
to -0.67 for distances above 800 Mpc (0.4 kpc/spaxel). This dis- 
tance effect would have implications in the derivation of the in- 
trinsic extinction, and subsequent properties, such as SFR, Ssfr, 
and the KS-law, in high-z star- forming galaxies, even in AO- 
based spectroscopy. If local LIRGs are analogues of the main- 
sequence star- forming galaxies at cosmological distances, the 
extinction values (Ay) derived from the observed emission lines 
in these high-z sources would need to be increased by a factor 
1.4 on average. 
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Fig. la! NGC 2369. Top panels show the SINFONI K band continuum emission, the observed maps (not corrected from extinction) of the lines Bry at 2.166//m 
and Br^ at 1.945 //m, together with the Ay map. The white contour englobes those spaxels above S/N = 4 considered to build the Ay map and distribution. The 
brightest spaxel of the K band continuum is marlced with a plus symbol (+), and has been used as reference to obtain the radial profile in the bottom right panel. The 
secondary nucleus, if present, is marked with a cross (x). Bottom left panel shows the the Av distribution of all valid spaxels (blue histogram) and the distributions 
of the those spaxels above (red) and below (yellow) the median value of the Br(5 S/N distribution. The relationship between the surface brightness of the lines and 
the Ay is shown in the central panel, where the points with Av<0 are outlined with a black contour Finally, the radial distribution of the extinction is shown in the 
bottom right panel, where the blue line represents the mean value of Ay for different radial bins and its error. The bins are obtained as the 1/30 of the total radial 
coverage of the map. 
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Fig. lb: NGC 3110. Same as Fig. la but for NGC 3110 
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Fig. Ic: NGC 3256. Same as Fig. la but for NGC 3256. Please note that the central spaxel lays outside the FoV since the nucleus 
was not observed in K-band. See text and Paper I for further details. 
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Fig. Id: ESQ 320-G030. Same as Fig. la but for ESQ 320-G030 
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Fig. le: IRASF 12115-4656. Same as Fig. labut for IRASF 12115-4656 
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Fig. If: NGC 5135. Same as Fig. la but for NGC 5135 
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Fig. Ig: IRASF 17138-1017. Same as Fig. labut for IRASF 17138-1017 
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Fig. Ih: IC 4687. Same as Fig. la but for IC 4687 
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Fig. li: NGC 7130. Same as Fig. la but for NGC 7130 
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Fig. Ij: IC 5179. Same as Fig. la but for IC 5179 
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Fig. 2a: IRAS 06206-6315. Top panels show the SINFONI K band continuum emission, the observed maps (not corrected from extinction) of the lines Paof at 
1.876//m and Bry at 2.166//m, together with the Ay map. The white contour englobes those spaxels above S/N = 4 considered to build the Ay map and distribution. 
The brightest spaxel of the K band continuum is marked with a plus symbol (+), and has been used as reference to obtain the radial profile in the bottom right panel. 
The secondary nucleus, if present, is marked with a cross (x). Bottom left panel shows the Av distribution of all valid spaxels (blue histogram) and the distributions 
of the those spaxels above (red) and below (yellow) the median value of the Bry S/N distribution. The relationship between the surface brightness of the lines and 
the Ay is shown in the central panel, where the points with Ay<0 are outlined with a black contour. Finally, the radial distribution of the extinction is shown in the 
bottom right panel, where the blue line represents the mean value of Ay for different radial bins and its error. The bins are obtained as the 1/30 of the total radial 
coverage of the map. For those objects with multiple components, the top inset shows the Ay radial profile of the system taking the main nucleus as the centre, 
whereas the bottom subpanel shows the radial profile obtained by extracting the profiles of each component separately and plotting them in the same reference frame. 
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Fig. 2b: IRAS 12112+0305. Same as Fig. 2a but for IRAS 12112+0305 
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Fig. 2c: IRAS 14348-1447. Same as Fig. 2a but for IRAS 14348-1447. 
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Fig. 2d: IRAS 17208-0014. Same as Fig. 2a but for IRAS 17208-0014. 
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Fig. 2e: IRAS 21130-4446. Same as Fig. 2a but for IRAS 21130-4446. 
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Fig. 2f: IRAS 22491-1808. Same as Fig. 2a but for IRAS 22491-1808. 
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Fig. 2g: IRAS 23128-5919. Same as Fig. 2a but for IRAS 23128-5919. 
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